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Stellar population gradients in early-type cluster galaxies 
with VIMOS IFU* 



oo 
O 

o 

(N 



(N 



6 



> 

in 

^J- 

en 

o 
oo 
o 



X 



T. D. Rawlef, Russell J. Smith, J. R. Lucey, A. M. Swinbank 

Department of Physics, Durham University, Durham DH1 3LE, United Kingdom 



2008 July 21 



ABSTRACT 

We present results from a pilot study of radial stellar population trends in early- 
type galaxies using the VLT VIMOS integral field unit (IFU). We observe twelve 
galaxies in the cluster Abell 3389 (z rs 0.027). For each galaxy, we measure 22 line- 
strength indices in multiple radial bins out to at least the effective radius. We derive 
stellar population parameters using a grid inversion technique, and calculate the radial 
gradients in age, metallcity and a-abundance. Generally, the galaxies in our sample 
have flat radial trends in age and [a/Fe], but negative gradients in [Z/H] (-0.20 ± 
0.05 dex). Combining our targets with two similar, long-slit studies to increase sample 
size, we find that the gradients are not correlated with the central velocity dispersion 
or iiT-band luminosity (both proxies for galaxy mass). However, we find that the age 
and metallicity gradients are both anti-correlated with their respective central values 
(to > 4er), such that galaxies with young cores have steeper positive age gradients, 
and those with metal-rich centres have strong negative [Z/H] gradients. 

Key words: galaxies: elliptical and lenticular, cD - galaxies: stellar content 



1 INTRODUCTION 

Stellar population gradients offer a valuable probe of the 
star formation and chemical evolution of early-typ e galax- 
ies. S i mple models of classic monolithic collapse (JLarsonl 
1 1974 ICarlberd I1984T ) can be ruled out by the rela- 
tive weakness of the observed metallicity gradients (e.g. 
iTamura fc Ohtal 120041 ). Generally, age gradients are ex- 
pected to be small for th e short collaps e times required 
for a-enhancement (e.g. iThomasI 1 19991 ), while the a- 
enhancement itself could show positiv e or negative gradients 
(|Pipino. D'Ercole fc Matteuccill2007f) . Within the hierarchi- 



19901 : ICarollo. Danziger fc Busonl Il993l; iMehlert et~ai] 
2003; Sanchcz-Blazqucz. Gorgas fc Cardiell 2006; 

Norris, Sharpies & Kuntschnerj 120061 : iReda et al.1 12007 



cal ACDM models (e.g . ICole et all l2000l : |Pe Lucia et~ai] 



120061 : iBower et al.l I2006D , merger gas content and feedback 
prescriptions can be investigated IjKawata fc Gibson! 120031 : 
iKobavashil 12004 ): mergers wou ld tend to dilute or erase 
Z/H gradients fe.g. lWhitelll980l ). although dissipation could 
trigger star formation at the core, creating age gr adients 
JBarnes fc Hernquistlll99ll : iMihos fc HernauistJll994 ) . 

Traditionally, observational studies of radial 
trends in early-type galaxies have employed long- 
slit spectroscopy in local cluster and field environ- 
ments (e.g. iGorgas. Efstathiou fc Aragon Salamancal 
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Sanchez-Blazquez et al.ll2007l ). In general, most studies find 
flat age gradients, small a-abundance gradients (either 
positive or negative), and large negative radial trends in 
metallicity. 

Long-slit spectrographs generally have a good optical 
throughput. However, the spatial data is immediately re- 
duced to one-dimension and is dependent on the slit orien- 
tation, with most of the light from the galaxy undetected. 
The use of an integral field unit (IFU), which enables both 
spectral and spatial coverage, allows the analysis of the en- 
tire galaxy. With careful treatment, light from the entire 
circumference can be binned to give radial trends unbiased 
by a slit orientation. 

One notable use of an IFU in this field is the Spec- 
trographic Areal Unit for R esearch on Optical Nebu- 
lae ( SAURON) project IJBacon et alj|200ll : Ide Zeeuw et ail 
2002). The instrument has good spatial resolution allow- 
ing the creation of maps for each measured parameter, 
particularly in kinematics. However, the spectral range of 
SAURON is limited to four useful line indices (H/3, Fe5015, 
Mgb5177, Fe5270 s ). Fro m spatial maps of the iron lines, 
iKuntschner et al.l i|2006h deduce that the metallicity gra- 
dients are generally negative and are roughly aligned with 
the light profiles. The H/3 maps are typically flat, with a 
strong central peak in a few galaxies, interpreted as recent 
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star formation activity in the core. Stellar population maps 
for the SAURON ea rly-type spiral sample are analysed in 
IPeletier etafl <|2007h . 

Here, we assess the feasibility of using the VLT Visible 
MultiObject Spectrograph (VIMOS) IFU for studying radial 
trends of stellar populations by observing a small sample of 
early-type galaxies in the nearby cluster Abell 3389 (z ~ 
0.027). With a large future study in mind, we develop the 
observing strategy and analysis techniques, comparing our 
results with those of previous work. The paper is organised 
as follows: Section [2] describes the sample, observations and 
data reduction; Section [3] details the measurement of the 
gradients in the absorption lines and in the stellar popu- 
lations; Section [4] analyses the stellar population gradients; 
Section [5] highlights our main conclusions. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 The sample and observations 

We observed a sample of 12 early- type galaxies from the 
nearby cluster Abell 3389 (cz « 8100 kms" 1 ). This study 
was designed as a feasibility assessment for the methodology 
and as such, the target cluster was selected based on its 
low-demand RA (~ 6 h ). Each galaxy is a confirmed cluster 
member and is morphologically classified as an early-type by 
iDresslerl Q1980, D80). A thumbnail image from the Digitized 
Sky Survejij is displayed in Figure [T] for each target, while 
Table [T] lists various fundamental parameters. 

The sample spans a range of three magnitudes in lu- 
minosity, and the galaxies have effective radii of 2.1 - 12.5 
arcsec (1.1 - 6.7 kpc). The effective radii (r c ) of the galaxies 
were measured from 2MASS J-band images. Due to the rel- 
atively small size of the target galaxies with respect to the 
2MASS PSF (FWHM ~ 3 arcsec), the apparent half-light 
radii required a correction. For this, we adopted the differ- 
ence be tween the half-lig ht radii of Sersic model fits (using 
GALFlT; |Peng et al]|2002h before and after PSF convolution. 
Our smallest object (D61) has an apparent half-light radius 
of 2.9 arcsec, which was corrected to 2.1 arcsec. The un- 
certainty in this correction, derived from the scatter in cor- 
rections for similar sized galaxies, is ±0.1 arcsec. For larger 
galaxies, the PSF correction is smaller, and we estimate that 
the listed effective radii have an uncertainty of better than 
10 per cent. 

Ten of the twelve target galaxies were previously ob- 
served through 2 arcsec diameter fibres as part of the 
NOAO Fundamental P lane Survey (NFPS; ISmith et all 
l2004lNelan et alj|2005r ). From the NFPS spectra, the only 
target to display line emission is D52 (EW[0 ill A5007] = 
0.55 A). D52 is a bright galaxy with a 'D'-type morphology, 
so the emission is most likely LINER-like rather than due to 
star formation. The emission line ratio [Oiii]/H/3 is greater 
than 3 for LINER-like objects and Hy/H/3 ~ 0.6, as inferred 
from nebular emission spectra. Therefore, in our worst-case 
galaxy, emission decreases the H/3 and H7F absorption in- 
dices by no more than 0.2 and 0.1 A respectively. Using the 
iThomas et al.l i|2003l . |2004 ) stellar population model grids, 
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Table 2. VIMOS instrument configuration. 



Mode 

Grism 

Field of view 

Spatial resolution 

Wavelength range 

Spectral resolution 

Dispersion 



IFU 

HR blue 

27 X 27 arcsec 

0.67 arcsec/fibre 

4150 - 6200 A 

2.1 A FWHM 

0.51 A/pixel 



we find that this causes an over-prediction in the ages of < 
5 per cent (< 1 Gyr for D52). 

The target galaxies were observed in October and 
November 2006 with the Visi ble MultiObject Spectrograph 
CVIMOS: ILeFevre et alJlioolh on the ESO VLT Melipal at 
Paranal. VIMOS was used in the Integral Field Unit (IFU) 
mode with the high resolution, blue grism. Table [5] de- 
tails the instrument configuration used. Observations were 
made during dark time, with an average seeing of 0.7 arc- 
sec (FWHM). Each galaxy was observed for 1730 seconds, 
achieving an average signal-to-noise ratio (SNR) ~ 18 A -1 
at 5000 A, for a 1 arcsec annulus at the effective radius 
(r c ). An offset sky frame with an exposure of 120 s was 
obtained for each target. Including overheads, a 1 hour ob- 
serving block was devoted to each galaxy. 



2.2 Data reduction and radial binning 

Initial data reduction, from raw detector output to sky- 
subtracted data-cube, was achieved using the ESO common 
pipeline (EsoRex) and the iraf routine imcombine. Each 
observation consists of three dithered exposures split into 
four quadrants, with an additional single exposure of four 
quadrants for the sky offset frame. For each quadrant, a sky 
background spectrum was calculated by taking the median 
over all spatial pixels in the offset exposure and then nor- 
malising, pixel-by-pixel, to the 5577 A line of the science 
frame. These sky spectra were then subtracted from the ap- 
propriate quadrants before performing a median combine 
of the three exposures. The resulting data-cubes were flux 
calibrated using the observed standard stars. 

For each galaxy, an ellipse was fit to the broadband, 
wavelength-collapsed ('white') image. This ellipse formed 
the basis for the pseudo-isophotic annuljj. Working from the 
centre of the galaxy outwards (in steps of 0.5 arcsec), indi- 
vidual fibre spectra were summed in elliptical annuli, with 
the bin boundaries based on two criteria: 

(i) annulus SNR > 18 A" 1 

(ii) annulus width (along major axis) Js 1 arcsec 

Two exceptions were adopted: the central bin is allowed 
a 0.5 arcsec semi-major axis provided that the minimum 
SNR is achieved; the outermost bin may have an SNR < 
18 A -1 if the radius is limited by the field of view. The 
targets have between three and seven radial bins, depending 
on effective radius and surface brightness. As an illustration, 



2 Named 'pseudo-isophotic' to stress that the annuli do not follow 
the true isophotes, but are based on the ellipse fit to the global 
galaxy shape. 
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Figure 1. Thumbnail images from the Digitized Sky Survey, showing a 100 arcsec square field around the target galaxies. The boxes 
correspond to the 27x27 arcsec VIMOS field of view, underlining how well matched our chosen galaxies are to the instrument. 



Table 1. Parameters of the target galaxies in Abell 3389. ID and morphology from Dressier {1980, D80). Positions from 2MASS. i?,-band 
magnitude from NED; /C-band magnitude from 2MASS. Effective circular radius (r e ) calculated from 2MASS J-band images. Central 
velocity dispersion, a (km s _1 ), as measured in Section 12.31 using the r c /3 aperture. 



ID 

(D80) 



NGC RA 

# (J2000) 



Dec 
(J2000) 



m R 
mag 



m K 
mag 



r e Morph. 

arcsec (D80) 



(km s 



D40 




06 21 59.7 


-64 59 19 


14.2 


11.16 


1.2 


SO 


178.8 ± 6.1 


D41 


2233 


06 21 40.1 


-65 02 00 


13.5 


10.47 


4.6 


SO 


184.6 ± 2.4 


D42 


2230 


06 21 27.5 


-64 59 34 


12.6 


9.66 


10.5 


E 


253.9 ± 3.1 


D43 




06 21 13.9 


-65 00 58 


14.6 


11.58 


3.1 


SO 


149.8 ± 4.4 


D44 




06 20 54.4 


-65 01 08 


14.0 


11.13 


3.8 


SO 


221.5 ± 4.5 


D48 




06 23 49.0 


-64 57 15 


14.8 


11.94 


2.9 


E 


144.6 ± 2.3 


D49 




06 23 07.4 


-64 55 50 


14.7 


11.64 


2.3 


E 


180.9 ± 3.9 


D52 


2235 


06 22 22.1 


-64 56 03 


12.4 


9.35 


12.3 


D 


252.1 ± 2.4 


D53 




06 22 04.8 


-64 57 36 


15.1 


12.18 


2.8 


E 


113.5 ± 4.4 


D60 




06 22 19.8 


-64 14 04 


13.5 


10.58 


1.5 


E 


194.0 ± 1.9 


D61 




06 21 58.3 


-64 53 44 


15.3 


12.37 


2.1 


SO 


146.6 ± 6.7 


D89 


2228 


06 21 15.6 


-64 27 32 


12.8 


9.79 


9.3 


SB0 


247.7 ± 2.8 
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Figure 2. Left panel: example of elliptical bins, showing broadband, wavelength-collapsed image from VIMOS for D60 with the adopted 
boundaries. Thick ellipse indicates the effective radius (r e ). Inset: total SNR within each annulus; r c shown as a dashed line. Right: 
spectra in the observed frame, with arbitrary offsets, for the five annuli within r c ; innermost annulus at the top. 
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Table 3. Quantitative comparison of VIMOS data to the previous 
NFPS line strength measurements. For offset comparisons, CNi, 
CN2, Mgi and Mg2 have been converted to A (from mags). 



Figure 3. Radial velocity (V; left panel) and velocity dispersion 
(ex; right panel) maps for the centre of D60 (effective radius indi- 
cated by the black ellipse). Both kinematic quantities in km s _1 , 
with the scale shown to the right of each panel. 



Figure [2] (left panel) shows the adopted bins for D60, along- 
side a histogram of the resulting SNRs (inset). For the same 
galaxy, total spectra from the five annuli within r c are given 
in the right panel of Figure [2] Additionally, we computed 
the spectra within three circular apertures: 

(i) 1 arcsec radius (for direct comparison to NFPS) 

(ii) effective radius, r c 

(iii) r c /3 (for galaxy-to-galaxy comparison) 

2.3 Velocity dispersion 

The stellar kinematics (V, a) of each galaxy were 
computed from the spectra of ea ch bin, using PPXF 
jCappellari fc Emsellemll2004l ) and the lVazdekisl l| 19991 ) evo- 
lutionary stellar population templates. These templates 
cover a similar wavelength range to VIMOS (we use ~4900- 
5500 A for fitting) and importantly, have a higher spectral 
sampling (FWHM = 1.8 A) than VIMOS HR blue. Errors 
were calculated using 51 Monte Carlo simulations, with the 
la intervals taken from the 8th and 42nd values after nu- 
merically sorting. 

Comparing velocity dispersion measurements from the 
1 arcsec radius circular aperture to NFPS, we find a good 
agreement with a median offset of 0.023 ± 0.008 dex (new 
values generally larger), and a scatter around the offset of 
0.019 dex. The scale also shows reasonable consistency, with 
a slope of 1.02 ± 0.06. 

For the galaxy with the highest total SNR within r c 
(D60), we created maps of the kinematics (V and a). Opti- 
mising for both SNR and spatial resolution, the data were 
re-binned into 16 segments of annuli, ranging from SNR~20- 
35 in the centre, to SNR~10-15 at the effective radius. The 
maps for V and a (Figure [3]) show that D60 has major 
axis rotation of ~110 kms -1 . Unfortunately, the total SNRs 
within r c for our other targets were insufficient for a simi- 
lar analysis. A one hour exposure with VIMOS is generally 
unsuited to the spatial mapping of galaxies at our target red- 
shift (previous studies of 2D absorption line kinematics, e.g. 
SAURON, observed galaxies an order of magnitude closer). 
The remainder of this paper adopts the radial binning de- 
scribed in Section \2. 21 



2.4 Line strength measurements 

Absorption-line strengths were measured from the flux- 
calibrated spectra at the VIMOS instrument resolution 





offset (A) 


X 2 1 DoF 


P(X 2 ) 




(1) 


(2) 


(3) 


H<5A 


-0.33 ± 0.40 


9.37 / 9 


0.60 


H<5F 


0.22 ± 0.22 


7.98 / 9 


0.46 


CNi 


-0.10 ± 0.30 


6.17 / 9 


0.28 


CN 2 


-0.04 ± 0.29 


4.77 / 9 


0.15 


Ca4227 


0.14 ± 0.20 


19.73 / 9 


0.98 


G4300 


0.12 ± 0.18 


11.31 / 9 


0.75 


H7A 


0.32 ± 0.36 


15.43 / 9 


0.92 


H7F 


0.03 ± 0.15 


13.62 / 9 


0.86 


Fe4383 


0.42 ± 0.39 


15.07 / 9 


0.91 


Ca4455 


0.08 ± 0.07 


4.28 / 9 


0.11 


Fe4531 


-0.06 ± 0.20 


20.83 / 9 


0.99 


Fe4668 


0.04 ± 0.17 


2.97 / 9 


0.03 


H/3 


0.13 ± 0.10 


15.33 / 9 


0.92 


Fe5015 


0.48 ± 0.30 


27.04 / 9 


1.00 


Mgi 


-0.12 ± 0.22 


7.68 / 9 


0.43 


Mg 2 


0.33 ± 0.17 


6.50 / 9 


0.31 


Mgb5177 


0.22 ± 0.10 


16.31 / 9 


0.94 


Fe5270 


0.37 ± 0.14 


24.07 / 9 


1.00 


Fe5335 


0.25 ± 0.10 


8.88 / 9 


0.55 


Fe5406 


0.26 ± 0.17 


22.84 / 6 


1.00 


Fe5709 


0.00 ± 0.12 


6.16 / 3 


0.90 


Fe5782 


0.04 ± 0.06 


9.71 / 9 


0.63 



NOTES - (1) Systematic offset of VIMOS measurements from 
NFPS. (2) x 2 ) given the systematic offset, and degrees of 
freedom (DoF). (3) Probability that x 2 is smaller than 
calculated. 



(~2.1 A FWHM; ~ 50 km s" 1 ) using INDEXI0. Uncertainties 
were calcul ated by indexf from the propagation of random 
errors (see ICardiel et al.lll998l ). We measured twenty- two 
atomic and molecular Lick indices in the VIMOS wavelength 
range (see first column of Table [3} . 

Model comparisons are performed in a relative sense, 
so we do not calibrate to the Lick instrumental re- 
sponse. However, the line indices do require corrections 
to account for velocity-broadening and instrumental res- 
olution. Here, we adop t the correction method used in 
ISmith. Lucev fc Hudson! i|2007l ). whereby the correction to 
zero velocity dispersion and to Lick resolution are performed 
in one step. Indices measured from a range of stellar popu- 
lation templates, smoothed to the Lick resolution, are com- 
pared to those measured from the same templates smoothed 
to the instrument resolution and further degraded to mimic 
various levels of velocity broadening. For each index, at a 
given velocity dispersion, this comparison gives a near lin- 
ear relation between the measured and corrected index. This 
technique removes the need to artificially broaden the ob- 
served galaxy spectra prior to index measurement, thereby 
preserving the non-uniform noise properties. 

Table [3] compares the central index measurements from 
VIMOS with those from NFPS for the 10 galaxies in com- 
mon. The VIMOS measurements are larger in most cases, 
with an average offset of 0.13 A. The x 2 > after accounting 



http://www.ucm.es/info/Astrof/software/indexf/ 
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Figure 4. Plots showing VIMOS 1 arcsec radius central region 
absorption line me asurements compared with those from NFPS 
IINelan et al.l 120051) . Solid line shows median offset; dashed line 
is equality. Values above each plot show the offset, \ 2 I DoF 
and P(x 2 ) respectively, as in Tabic [3] Given the actual scatter, 
the formal errors on H/3 and Mgb5f77 appear under-estimated, 
whereas for Fe4668 they appear over-estimated. 
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Figure 5. Line-strength index gradients versus central logo - . The 
solid grey line gives the median gradient, which is also displayed 
at the top of each panel. The dashed black line indicates zero 
gradient. 



for the sample offsets, suggests a reasonable agreement, al- 
though the errors on the indices are, in some cases, over- or 
under- estimated by up to 20 per cent (see Figure [4}. 

In AppendixfA]we tabulate aperture corrections derived 
from the integrated line strengths within two circular aper- 
tures (with radii of r c and r e /3). 



3 MEASUREMENT OF GRADIENTS 

Gradients (for line strengths and stellar population param- 
eters), along with their errors, are calculated using an un- 
weighted linear least-squares fit to all data points within the 
effective radius of the galaxy, with radii in log(r/r c ). Our 
choice of using log-space is not optimal for every galaxy as 
some targets would certainly be better fit if it were linear, 
i.e. r/r a . The necessity of consistency dictates the use of the 
same radial units in every case, and the majority of the tar- 
gets/parameters prefer a fit in log-space. Also in favour of 
log(r/r e ) is its adoption by previous studies of radial trends. 
The un- weighted fit ensures sensible results for the minority 
of cases which would be better suited to linear radial units. 

Three of the four smallest galaxies (D48, D53, D61; ef- 
fective radius < 3 arcsec) have insufficient radial bins with 
adequate SNR for determining gradients, and as such have 
not been included in the analysis. The remaining 'small' 
galaxy (D49) has a higher surface brightness, allowing three 
annuli with SNR ^ 18 within the effective radius. 

We first discuss the gradients of the absorption indices 
(Section 13. ip . In Section 13.21 we calculate the median line 
strength gradients, from which we infer radial trends in the 
stellar population parameters, averaged over all galaxies. Fi- 



nally, a grid inversion technique is used to derive the stel- 
lar population gradients for each individual galaxy (Section 

nnj. 



3.1 Line-strength gradients 

The line-strength gradients of the 22 measured indices, for 
each galaxy in our sample, are tabulated in Appendix iBl 

As an illustration, Figure [5] shows the measured line- 
strength gradients (against velocity dispersion) for eight of 
the indices. The median gradient for the lines is given in each 
panel. Interestingly, line indices which to first order follow 
the same physical property, have quite different gradients. 
For example, the two 'age-tracing' Balmer lines show very 
different gradients: H7F generally has a positive gradient 
(0.50 ± 0.24) while H/3 has a flat profile (-0.04 ± 0.24). 
This can be explained by a higher sensitivity to metallicity 
of the H7F index. 



3.2 Average stellar population gradients 

By exploring the response of each index to age, metallicity 
[Z/H] and a-abundance [a/Fe], we can describe the average 
radial stellar population trends in our sample. 

The linear response of each index to the underly- 
ing stellar pop u lation param eters were calculated from the 
iThomas et ail |2003l . J2004 ) model grids. The responses 
quantify the change in each index for a decade c hange in each 
of age , Z/H and a/Fe, derived as described in ISmith et al.l 
\200Jt ). 
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Figure 6. H/3 gradient versus H7F gradient. The solid line gives 
the predicted index gradients for a simple model with a [Z/H] 
gradient only (no radial trend in age or [a/Fe]). The two main 
outliers, presumably the galaxies with strongest age gradients, are 
labelled. 



Table 4. Mean stellar population gradients. The first two rows 
show gradients for our VIMOS sample, as derived by two meth- 
ods: (1) from the median line strength gradients and index re- 
sponses (Sec. 13.21 1: (2) individual galaxy grid inversion s (Sec . 
331 Rows 3 and 4 give gradien ts from iMehlert et"al] (120031) 
and lSanchez-Blazquez et al.l (120071 ) respectively, and the final row 
shows our combined sample (see Sec. l4.lt . 





log(age/Gyr) 


[Z/H] 


[a/Fe] 


This study (1) 
This study (2) 


0.07 ± 0.06 
0.08 ± 0.08 


-0.21 ± 0.05 
-0.20 ± 0.05 


0.06 ± 0.03 
-0.02 ± 0.03 


M03 

SB07 

Combined 


0.03 ± 0.04 
0.15 ± 0.10 
0.06 ± 0.03 


-0.12 ± 0.02 
-0.31 ± 0.04 
-0.18 ± 0.02 


-0.01 ± 0.01 

0.01 ± 0.02 

-0.01 ± 0.01 



ity gradients, but fiat radial trends for log(age/Gyr) and 
a-abundance. 



Figure [6] compares the H7F and H/3 Balmer line gra- 
dients for each galaxy. Given the derived responses of the 
Balmer lines, the relation for a metallicity gradient only (no 
radial trend in age or [a/Fe]), is shown. Two galaxies (D43, 
D49) have a large offset from the predicted Balmer gradi- 
ents of this simple model. The H/3 response to [a/Fe] is too 
small to explain the discrepancy, requiring a very large neg- 
ative gradient. This suggests that these galaxies have the 
strongest age gradients, and given the negative age response 
of the Balmer lines, the age gradient is expected to be pos- 
itive. We confirm this interpretation in Section [4] 

We consider several triplet combinations of indices, each 
including a Balmer line, an iron line and Mgb5177. From the 
mean gradients of the three indices, we use the responses 
to solve for the average gradients in age, metallicity and 
a-abundance (Figure [7|). The left panel indicates that the 
derived age and metallicity gradients are dependent on the 
Balmer line/iron line combination, with a scatter in the di- 
rection of the usual age-met allicity degeneracy. The right 
hand panel highlights that the [a/Fe] gradients tend to com- 
pensate for the age in these models, but also shows that the 
[a/Fe] gradient is generally independent of which Balmer 
line is used. The triplets including Fe4383 and Fe4531 have 
the largest uncertainties, due to the large errors for the me- 
dian line strength gradients used in the stellar population 
reconstruction. 

The dependence of the derived metallicity gradi- 
ents on the indices used w as previously discussed in 
ISanchez-Blazquez et al.l (|2006l ). These authors found that 
the index pair Fe4668-H/3 gave a metallicity gradient ~0.4 
dex steeper than Fe4383-H/3, a significantly larger difference 
in metallicity gradient than seen in our equivalent triplets 
(the index pair plus Mgb5177), where the difference is ~0.07 
de x. This underlin es that the inclusion of a-enhancement via 
the lThomas et al.l models and Mgb5177 yields a higher level 
of consistency. 

The mean stellar population gradients for the sample, 
as estimated from the triplets of lines, are listed in the first 
row of Table [4] Generally, they support negative metallic- 



3.3 Individual stellar population gradients 

We derive stellar population parameters for each ra- 
dial bin of each galaxy using a grid inversion technique 
IjProctor. Forbes fc Beaslevl |2004| ; Smith et al. in prep.). 
This involves a non-linear % 2 minimisation over several in- 
dices, using a cubic-spline interpolation in age-[Z/H]-[a/Fe], 
and a linear extrapolation for points lying outside the model 
grid. 

We use the lThomas et al.l (|2003l . 120041 ) models together 
with eight measured indices: two Balmer lines (H7F, H/3), 
five iron lines (Fe4384, Fe4531, Fe4668, Fe5270, Fe5335) and 
the a/Fe sensitive Mgb5177. The shorter wavelength Balmer 
indices (H<5A, H<5F) are excluded for having low SNR due 
to the instrumental response, and H7A is redundant as it 
mostly duplicates the response of H7F. The iron lines at the 
red end (Fe5406, Fe5709, Fe5782) are affected by the 5577 A 
sky line and other sky features at the redshift of our cluster. 

Stellar population gradients for each galaxy are calcu- 
lated in an analogous manner to the line strength gradients 
in Section r3.ll The stellar population gradients are listed in 
Appendix Table IC1I and plotted for individual galaxies in 
Figures ICi1-IC3l 



4 RESULTS AND DISCUSSION 

The distributions of log(age/Gyr), [Z/H] and [a/Fe] gradi- 
ents for our sample of A3389 galaxies are shown in the top 
row of Figure [8] The mean age and a-abundance gradients 
are consistent with flat radial trends (0.08 ± 0.08 and -0.02 
± 0.03 respectively), while the mean metallicity gradient is 
-0.20 ± 0.05. The age and metallicity gradients are in agree- 
ment with those derived from the mean index gradients (in 
Sec. I3.2[l . as summarised in Table [4] The a-abundance gra- 
dients calculated by the two methods marginally disagree. 
The two galaxies in the most positive age gradient bin (see 
top row of Fig. [8} are D43 and D49, which supports the 
prediction from the Balmer line plot in Fig. [6] 
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Figure 7. Average stellar population gradients over the nine galaxies, derived from mean gradients of various triplets of line indices. 
All triplets incorporate a Balmer line, an iron line and Mgb5177. Points are represented by the iron line wavelength. Black labels indicate 
H7F; grey indicates H/3. The sensitivity of the stellar population gradients to the choice of indices is evident. 
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Figure 8. Distribution of (left-to-right) log(age/Gyr), [Z/H] and 
[a/Fe], Shown from top to bottom are our VIMOS data, M03, 
SB07 and the combined sample (see text for details). Mean gra- 
dients are marked with a dotted line and displayed at the top left 
of each panel. 



authors rebinned the long-slit spectra to achieve a SNR of 
at least 30 out to the effective radius. Stellar population pa- 
rameters were derived by an index-pair iteration technique, 
using H/3, Mgb51 77 and (Fe) = (Fe5270 +Fe5335)/2, in con- 
junction with the lThomas et al.l (J2003J) models. SB07 anal- 
ysed a sample of 11 local early- type galaxies from a mix of 
environments (field, group, Virgo cluster). The authors re- 
binned the long-slit spectra along the axis to obtain a mini- 
mum SNR ~ 50 out to 2r c , and employed the lThomas et al.l 
models, using a % 2 minimisation technique on 19 Lick in- 
dices. The second and third rows of Figure [8] show the gra- 
dient distributions from these studies. 

Our results agree well with those from the previous 
studies (see Table [3J. In the combined sample, the in- 
dividual points are not weighted by errors, so the to- 
tals are dominated by the larger (but noisy) M03 study. 
We have not included t he ea rly-type galaxy sample from 
ISanchez-B fazqucz ct al. (2006) in our comparisons as they 
use the lVazdekisI ( 19991 . and in prep.) stellar population mod- 
els, which do not explicitly include a-enhancement. 

We split the combined data from the three samples (VI- 
MOS, M03, SB07) into the two traditional early-type mor- 
phology classes (elliptical and SO). Using a Student's t-test 
to analyse the distributions, there is no significant difference 
between the classes in any of the stellar population param- 
eters. 



4.1 Comparison studies 

Stellar population gradients may be affected by numerous 
variables (environment, morphology etc.), which is likely 
to result in a large intrinsic scatter within the measured 
and derived parameters. For small samples, this can con- 
tribute to misleading conclusions. We therefore analyse our 
results alongside d ata from two previous , long-slit studies of 
early- type galaxies: iMehlert et al.l (120031 . hereafter M03) and 
ISanchez-Blazquez et al.l (|2007l . SB07), chosen due to their 
use of the same stellar population models. M03 observed a 
sample of 35 galaxies in the Coma cluster, all morphologi- 
cally typed as early-type by D80. For their analysis, these 



4.2 Stellar population gradients versus central 
log o 

In the following sections, we compare the gradients to central 
values. The central properties of our VIMOS galaxies are 
from circular apertures with a radius one third that of the 
galaxy effective radius. This allows more meaningful galaxy- 
to-galaxy comparisons than if a fixed radius circle were used 
in all cases. The long-slit comparison studies measure central 
values as follows: M03 uses an average along the major axis 
inside 0.1r e /\/l — e , where e is the ellipticity (from optical 
imaging); SB07 employs an aperture size of 1.5 arcsec x 
r./8. 

We explore the relation between stellar population gra- 
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Figure 9. Gradient of log(age/Gyr) (left) and [Z/H] (right) ver- 
sus central logo - . Rows for our VIMOS galaxies, M03, SB07 and 
the combined sample are shown. Bottom row: solid line is an un- 
weighted linear best fit to all points; the slope (and significance) 
is shown in the top right of the panel. The two galaxies with 
the largest positive age gradients are GMP3201 (from M03) and 
NGC2865 (SB07). 



dients and the central velocity dispersion (logo - , a proxy 
for galaxy mass). The a-enhancement of g alaxies suggest s 
the co-requisite of small age gradients (e.g. lThomaslll999f) . 
while the metallicity gra dients are dependent on merge rs and 
feedback processes (e.g. IWhitelll980l ; lKobavashill2004l ) . In a 
simple model, enrichment gradients depend on the galactic 
wind efficiency, which varies with the size of the potentia l 
well and hence the galaxy mass (|Arimoto fc Yoshiilll987h . 
Figure [5] shows the log(age) (left panel) and [Z/H] (right 
panel) gradients versus the central log a. There appears to 
be no correlation of age gradient with log a (slope: -0.32 ± 
0.30). The same is true for the [a/Fe] gradient (not shown), 
which has a slope of -0.04 ± 0.08 in the combined sample. 
The [Z/H] gradient weakly correlates with log a, having a 
slope of 0.26 ± 0.17 in the combined sample (although not 
not statistically significant at 1.5a). 

SB07 claim a discontinuity at logo — 2.3 (a = 200 
km s _1 ) in the trend with metallicity, which can be seen 
in the right-hand panel, third row of Fig. [9] They specu- 
late that physically, this could mark the transition between 
'wet' and 'dry' mergers, caused by the interplay between 
in-flowing gas and the feedback processes (see also iBinnevi 
l2004i : lFaber et alj|2007j ). Unfortunately, our VIMOS sample 
does not include galaxies of low enough mass to test this 
result. However, the turning point is not seen in the larger, 
but noisier M03 data (which also dominates the combined 
plot ; bottom row i n Fig . [5J. 

iForbes et al. (2005), using a sample from M03 and 



ISanchez-Blazquea (J2004), claim a 99 per cent probability 
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Figure 10. Metallicity gradient versus K band absolute magni- 
tude (2MASS) for the combined sample (VIMOS+SB07+M03). 
Solid line is an un-weighted linear best fit to all points; the slope 
(and significance) is shown in the top right of the panel. 



of a correlation be tween [Z/H] gradient and logo". Addition- 
ally, IForbes et al.l find a 99.6 per cent probability of a corre- 
lation between [Z/H] gradient and the A'-band luminosity (a 
proxy for stellar mass; from 2MASS), whereas in our com- 
bined sample we find no significant trend (Figure [lOj . 

4.3 Stellar population gradients versus central 
values 

In Figure [11] (left-hand panels) we explore the age gradient 
versus the central value of the age. The bottom panel, which 
combines the three data sets, shows the best fitting line to 
the data. The SB07 data illustrates why an un-weighted fit is 
used throughout this study: several of the galaxies from that 
study have quoted age gradients of 0.000 ± 0.000 as a con- 
sequence of fitting to discrete model data. In the combined 
sample, the age gradient has a significant correlation with 
the central value (4.6 a), with a calculated slope of -0.46 ± 
0.10. The un-weighted fit is dominated by the M03 sample, 
and the magnitude of the slope may be unduly affected by 
the two centrally young galaxies with large positive gradi- 
ents (GMP3201 from M03; NGC2865 from SB07). Without 
these galaxies, the slope (and significance) is reduced to - 
0.20 ± 0.08 (2.4 cr). 

The central panels of Figure [11] shows the analogous 
plot for metallicity. The un-weighted linear fit to the com- 
bined sample gives a significant trend (4.2 a) with a slope of 
-0.38 ± 0.09. This general trend is readily seen in each of 
the individual datasets, and is noted by SB07 (albeit at 2 cr 
significance in their small sample) . Using central values from 
an independent study, SB07 conclude that this trend is not 
due to the correlation of errors. The correlation in our com- 
bined sample would be strengthened (slope of -0.56 ± 0.12; 
4.8 a) by removing the three galaxies from M03 (GMP3329, 
GMP2921, GMP3561) with anomalously metal-rich centres; 
M03 notes that these galaxies lie outside the model grid 
leading to large uncertainties in extrapolation. 

Generally, no trends were found between the a- 
abundance gradient and central properties of the galaxy. 
The right-hand panels of Figure [TT1 shows the lack of a trend 
between central [a/Fe] and its gradient. SB07 tentatively in- 
voke the dichotomy at log a = 2.3 for central a-abundance 
versus the gradient (third row of right-hand panel; open cir- 
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Figure 11. Gradient versus central value for log(age/Gyr) (left-hand column), [Z/H] (central column) and [a/Fe] (right-hand column). 
Layout as in Fig. [9] Additionally, in the right-hand column, open circles for log a > 2.3 else closed circles. 



cles for logo - > 2.3, else closed circles). However, the di- 
chotomy is not recovered in our data, or in M03. The ap- 
parent division may be a consequence of the small sample 
size, with the well known correlation between central ve- 
locity dispersion a nd central [a/Fe] (e.g. iGraves et al.ll2007l; 
ISmith et a l. 2007) ensuring that on average, the high mass 
galaxies lie further to the right in the plot. 

In summary, the observed trends imply that galaxies 
with young cores have steeper positive age gradients, and 
those with metal-rich centres have strong negative [Z/H] 
gradients. Assuming stellar population gradients behave lin- 
early out to the effective radius and continue to do so be- 
yond it, the trends suggest that all galaxies have a similar 
age and metallicity (old and [Z/H] w 0) at ~ 2-2.5 r c . Of 
the 5 galaxies in SB07 which measure [Z/H] to 2 r c , four are 
approximately solar metallicity at this radius and the other 
is metal-poor. Verification of this result is beyond the scope 
of our current VIMOS data. 



4.4 Correlation of errors 

For any particular radial bin, the age and metallicity are 
not truly independent. The well-known degeneracy between 
age and [Z/H] is to an extent broken when fitting model 
spectra to multiple line strength indices. Here, however, the 
gradients in these parameters suffer from a correlation of 
errors. 

To investigate the importance of correlated errors, we 
perform 10 3 Monte Carlo simulations for each galaxy in our 



A3389 sample. The realizations consists of a simulated age 
and metallicity at each radial bin, using a bivariate nor- 
mal distribution derived from the error ellipse associated 
with the measured points. The gradients in each parameter 
are then calculated as with the measured data. Figure [12] 
(upper panel) shows the measured data as points and the 
simulated values represented by 1 a confidence ellipses. The 
distribution of the points is parallel to the error ellipses, 
which is reminiscen t of the age-mass-met alicity ('Z') plane 
l|Trager et al.ll2000l ). 



To test whether any of the scatter is intrinsic, the mea- 
sured data points are projected onto the line of best fit, 
using the simulated ellipse major axes as a quantification of 
the measurement error. Hence, the distribution is reduced 
to one-dimension along the diagonal (Fig. 1121 lower panel). 
Taking the null hypothesis of no intrinsic scatter, the % 2 is 
23.69 (8 degrees of freedom) with the probability of obtain- 
ing this distribution from the hypothesis calculated as Q(x 2 ) 
= 0.003. 



Although there is a sizeable correlation of errors, we 
conclude that there is also a real intrinsic scatter in the 
same direction. Hence, the correlation we find between the 
central value and gradient of age is not merely a reflection 
of the equivalent trend in metallicity, and vice versa. The 
trends between the central value and gradient in age and 
metallicity are both real. 
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Figure 12. Upper panel: measured age and mctallicity gradients 
shown as points, with the la confidence ellipses obtained from 10 3 
Monte Carlo simulations for each of our sample galaxies. Line of 
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of points projected onto the best fit line (shown with an arbitrary 
y-offsct for clarity). Error bars derived from the simulations. The 
X-squared / degrees of freedom and Q (x 2 ) are shown to the lower 
left. The presence of a real intrinsic scatter is evident. 



5 CONCLUSIONS 

Using the VIMOS integral field unit, we observed a sample 
of early-type galaxies in the local cluster Abell 3389. From 
the resulting data-cube, we measured 22 line indices in ra- 
dial bins out to the effective radius, and successfully derived 
stellar population gradients. Our key findings are that early- 
type galaxies have: 

• strong negative gradients in metallicity (~ -0.2 dex per 
decade in radius); 

• radial trends in age and a-abundance that are consis- 
tent with zero; 

• stellar population gradients not significantly correlated 
with the central velocity dispersion or A'-band luminosity 
(proxies for total galaxy and stellar mass respectively); 

• a strong anti-correlation between the central values and 
gradients of both age and [Z/H]. 

In a future paper we will report on an additional sam- 
ple of 19 early-type galaxies in the Shapley Supercluster, 
observed with the VIMOS IFU. 
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APPENDIX A: DERIVED APERTURE 
CORRECTIONS 

We use the integrated line strengths from two circular aper- 
tures, with radii of r e /3 and r c , to derive aperture correc- 
tions. These corrections should aid index line strength trans- 
lations between apertures of different sizes; for example, be- 
tween different multi-object spectrographs. 

For each index, the median corrections [I(r c ) — I(r e /3)] 
over the 9 VIMOS galaxies included in the gradient anal- 
ysis (i.e. excluding D48, D53, D61) are presented in Table 
IA1I Corrections are also tabulated for the velocity disper- 
sion (as calculated in Section ^. 3[) and the stellar population 
parameters (as derived in Section f3 . 3 f) . 



Table Al. Median index line strength aperture corrections for 
circular apertures, I(r e ) — I(r e /3), for the 9 VIMOS galaxies with 
gradient data (i.e. excluding D48, D53, D61). Values for line in- 
dices in A, except CNi, CN2, Mgi and Mg2 which are in mags. 
Ql &: Q3 give the first and third quartiles of the data. Also tab- 
ulated are the equivalent values for the velocity dispersion and 
stellar population parameters. 





median correction 


Ql 


Q3 


H<5A 


-0.03 ± 0.48 


-0.91 


0.91 


HoF 


-0.12 ± 0.34 


-0.38 


0.53 


CNi 


0.01 ± 0.01 


-0.03 


0.03 


CN 2 


0.01 ± 0.01 


-0.02 


0.02 


Ca4227 


0.01 ± 0.15 


-0.06 


0.10 


G4300 


-0.03 ± 0.35 


-0.20 


0.94 


H7A 


0.09 ± 0.33 


-0.53 


0.37 


H7F 


-0.32 ± 0.24 


-0.37 


0.30 


Fe4383 


0.37 ± 0.35 


-0.29 


1.01 


Ca4455 


-0.07 ± 0.11 


-0.31 


0.03 


Fe4531 


0.11 ± 0.16 


-0.34 


0.25 


Fe4668 


-0.53 ± 0.33 


-0.99 


-0.29 


H/3 


-0.03 ± 0.09 


-0.19 


0.08 


Fe5015 


0.04 ± 0.23 


-0.24 


0.45 


Mgi 


-0.01 ± 0.01 


-0.01 


-0.01 


Mg 2 


-0.01 ± 0.01 


-0.03 


-0.01 


Mgb5177 


-0.33 ± 0.11 


-0.51 


-0.06 


Fe5270 


-0.10 ± 0.19 


-0.21 


0.05 


Fe5335 


-0.12 ± 0.09 


-0.36 


-0.05 


Fe5406 


-0.12 ± 0.10 


-0.20 


0.21 


Fe5709 


-0.04 ± 0.29 


-0.14 


0.21 


Fe5782 


-0.07 ± 0.05 


-0.10 


0.12 


logo- 


-0.02 ± 0.02 


-0.04 


0.00 


log(age/Gyr) 


0.03 ± 0.06 


0.00 


0.13 


[Z/H] 


-0.14 ± 0.04 


-0.17 


-0.09 


[a/Fe] 


-0.04 ± 0.03 


-0.12 


0.00 



APPENDIX B: LINE STRENGTH GRADIENTS 



Table [BT1 details the measured kinematics and lines strength 
gradients for nine of our galaxies (excluding D48, D53, D61 
which have insufficient radial bins for gradient determina- 
tion; see Section [3J) . 



APPENDIX C: 
GRADIENTS 



STELLAR POPULATION 



Figures IC1I IC2I and IC3I show the stellar population param- 
eters versus radius, in log(r/r c ), to the effective radius for 
all target galaxies. Table [Cll lists the measured gradients for 
loger, log(age/Gyr), [Z/H] and [a/Fe], for the nine galaxies 
analysed. 
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Table Bl. Absorption index gradients for the galaxies excluding D48, D53, D61. Identifications from iDresslerl (J1980J). 



ID 


H5A 


H<5F 


CNi 


CN 2 


Ca4227 


G4300 




H 7 A 


H 7 F 


Fe4383 


Ca4455 


Fe4531 


Fe4668 




H/3 


Fe5015 


Mgi 


Mg 2 


Mgb5177 


Fe5270 




Fe5335 


Fe5406 


Fe5709 


Fe5782 






D40 


1.016 ± 0.595 


0.532 ± 0.629 


0.611 ± 0.930 


-0.276 ± 0.168 


0.048 ± 0.028 


0.700 ± 0.489 




-0.338 ± 0.754 


0.602 ± 0.401 


0.108 ± 0.060 


-0.522 ± 0.258 


-1.127 ± 0.853 


-3.230 ± 2.231 




0.697 ± 0.742 


0.448 ± 0.775 


1.034 ± 0.069 


0.955 ± 0.356 


0.235 ± 0.093 


-0.016 ± 0.020 




-0.466 ± 0.143 


-0.033 ± 0.729 


-0.161 ± 0.089 


0.151 ± 0.100 






D41 


-2.029 ± 0.977 


0.350 ± 0.495 


1.329 ± 0.744 


1.697 ± 0.804 


0.252 ± 0.106 


2.601 ± 0.382 




-1.646 ± 1.225 


-0.353 ± 0.341 


-1.235 ± 0.429 


-0.416 ± 0.159 


-1.182 ± 0.334 


-1.861 ± 0.921 




0.350 ± 0.441 


1.546 ± 0.670 


-1.154 ± 0.501 


-0.714 ± 0.168 


-0.123 ± 0.225 


-0.934 ± 0.245 




-0.232 ± 0.355 


0.183 ± 0.499 


-0.601 ± 0.215 


0.009 ± 0.176 






D42 


-1.565 ± 1.607 


-2.049 ± 1.331 


0.740 ± 0.887 


-0.985 ± 0.437 


1.153 ± 0.487 


0.768 ± 0.770 




0.618 ± 0.576 


1.379 ± 0.265 


1.312 ± 0.562 


0.184 ± 0.159 


0.308 ± 0.209 


-0.957 ± 0.440 




-0.039 ± 0.081 


-0.793 ± 0.218 


-1.008 ± 0.315 


-0.053 ± 0.191 


-0.352 ± 0.272 


0.105 ± 0.245 




-0.067 ± 0.135 


-0.597 ± 0.208 


-0.300 ± 0.180 


0.130 ± 0.176 






D43 


1.300 ± 0.780 


0.374 ± 0.133 


-0.739 ± 1.018 


-0.620 ± 1.265 


0.390 ± 0.376 


-1.345 ± 1.393 




2.456 ± 1.095 


0.504 ± 0.565 


-1.238 ± 0.501 


-0.203 ± 0.115 


-0.260 ± 0.079 


-1.129 ± 0.263 




-1.026 ± 0.107 


-1.040 ± 0.397 


-1.876 ± 1.084 


-0.033 ± 0.443 


-0.441 ± 0.146 


-0.470 ± 0.176 




-0.412 ± 0.515 


0.222 ± 0.658 


0.001 ± 0.377 


-0.459 ± 0.037 






D44 


-0.008 ± 0.500 


-0.510 ± 0.315 


0.162 ± 0.524 


-0.302 ± 0.313 


-0.341 ± 0.072 


0.347 ± 0.349 




-0.257 ± 0.348 


-0.255 ± 0.190 


-0.230 ± 1.348 


-0.442 ± 0.277 


0.638 ± 0.469 


0.226 ± 0.541 




-0.261 ± 0.111 


-0.337 ± 1.029 


-1.186 ± 0.084 


-1.055 ± 0.804 


-1.128 ± 0.342 


-0.472 ± 0.261 




-0.490 ± 0.248 


-0.472 ± 0.266 


0.244 ± 0.068 


-0.550 ± 0.308 






D49 


-2.312 ± 0.517 


-1.392 ± 0.008 


0.798 ± 0.613 


0.717 ± 0.805 


-0.133 ± 0.015 


-0.810 ± 0.667 




0.365 ± 0.598 


0.184 ± 0.302 


1.263 ± 1.230 


-0.451 ± 0.087 


-0.638 ± 0.182 


-1.424 ± 0.211 




-0.862 ± 0.091 


0.019 ± 0.884 


-0.606 ± 1.581 


-0.926 ± 1.189 


-0.186 ± 0.751 


0.281 ± 0.444 




-0.952 ± 0.109 


0.443 ± 0.180 


1.178 ± 0.222 


0.232 ± 0.087 






D52 


-0.685 ± 0.402 


0.435 ± 0.236 


0.245 ± 0.291 


0.674 ± 0.309 


0.285 ± 0.135 


-0.854 ± 0.553 




1.400 ± 0.401 


0.934 ± 0.144 


0.947 ± 0.391 


-0.769 ± 0.235 


-0.630 ± 0.194 


-0.130 ± 0.392 




0.419 ± 0.274 


-0.302 ± 0.278 


-1.734 ± 0.289 


-1.190 ± 0.294 


-0.380 ± 0.123 


-0.053 ± 0.130 




-0.256 ± 0.253 


-0.693 ± 0.067 


-1.524 ± 0.469 


-0.136 ± 0.052 






D60 


1.133 ± 0.606 


0.776 ± 0.393 


-2.291 ± 0.503 


-2.690 ± 0.414 


-0.586 ± 0.439 


-0.373 ± 0.338 




0.580 ± 0.302 


0.622 ± 0.185 


-0.275 ± 0.342 


-0.063 ± 0.109 


0.026 ± 0.200 


-2.185 ± 1.117 




-0.028 ± 0.074 


-0.494 ± 0.301 


-1.541 ± 0.224 


-2.228 ± 0.782 


-0.901 ± 0.289 


-0.443 ± 0.115 




-0.944 ± 0.306 


-0.331 ± 0.106 


-0.264 ± 0.149 


-0.518 ± 0.057 






D89 


-0.429 ± 0.672 


-1.220 ± 0.365 


0.686 ± 0.481 


0.679 ± 0.399 


0.561 ± 0.279 


-0.464 ± 0.465 




-1.588 ± 0.477 


-0.101 ± 0.359 


1.854 ± 0.792 


0.148 ± 0.151 


-1.762 ± 0.369 


-0.819 ± 1.151 




-0.222 ± 0.192 


1.681 ± 0.389 


-2.398 ± 0.545 


-1.716 ± 0.670 


-1.171 ± 0.436 


-0.588 ± 0.268 




-0.584 ± 0.096 


-0.040 ± 0.139 


0.206 ± 0.134 


0.158 ± 0.125 







Table CI. Kinematic and stellar population gradients for the galaxies excluding D48, D53, D61. Identifications from IDresslerl 1119801) . 



ID 


gr ad (logo - ) 


grad(log(age/Gyr)) 


grad([Z/H]) 


grad([a/Fe]) 


D40 


0.039 ± 0.018 


0.094 ± 0.159 


-0.241 ± 0.076 


0.096 ± 0.038 


D41 


-0.006 ± 0.063 


0.079 ± 0.155 


-0.168 ± 0.109 


0.075 ± 0.031 


D42 


-0.051 ± 0.009 


-0.253 ± 0.177 


0.022 ± 0.068 


-0.098 ± 0.040 


D43 


-0.059 ± 0.048 


0.466 ±0.116 


-0.354 ± 0.069 


0.004 ± 0.107 


D44 


-0.171 ± 0.010 


0.004 ± 0.110 


-0.167 ± 0.133 


-0.112 ± 0.091 


D49 


-0.054 ± 0.038 


0.425 ± 0.172 


-0.326 ± 0.004 


-0.005 ± 0.147 


D52 


0.037 ± 0.016 


-0.259 ± 0.099 


0.006 ± 0.074 


-0.044 ± 0.034 


D60 


-0.009 ± 0.032 


0.067 ± 0.199 


-0.369 ± 0.149 


0.011 ± 0.088 


D89 


-0.164 ± 0.037 


0.064 ± 0.083 


-0.239 ± 0.129 


-0.146 ± 0.091 
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Figure CI. log(age/Gyr) versus log(r/r e ). Grey dotted lines show the bin boundaries. Black dashed line is at the effective radius. The 
fit to point s (solid black line) within r c is un- weighted, with the gradient given in the top left of each panel. Name and morphological 
type (from iDresslerllf 9801 ) are given to the bottom right and left respectively. Gradients are not calculated for D48, D53, D61. 
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Figure C2. [Z/H] versus log(r/r c ). Layout as in Fig. lCll 
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Figure C3. [a/Fe] versus log(r/r c ). Layout as in Fig. lCll 



